Stability of conductance oscillations in monatomic sodium wires 
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We study the stability of conductance oscillations in monatomic sodium wires with respect to 
structural variations. The geometry, the electronic structure and the electronic potential of sodium 
wires suspended between two sodium electrodes are obtained from self-consistent density functional 
theory calculations. The conductance is calculated within the framework of the Landauer-Biitttiker 
formalism, using the mode- matching technique as formulated recently in a real-space finite-difference 
scheme [Phys. Rev. B 70, 195402 (2004)]. We find a regular even-odd conductance oscillation as 
a function of the wire length, where wires comprising an odd number of atoms have a conductance 
close to the quantum unit Go = e^ /nh, and even-numbered wires have a lower conductance. The 
conductance of odd-numbered wires is stable with respect to geometry changes in the wire or in 
the contacts between the wire and the electrodes; the conductance of even-numbered wires is more 
sensitive. Geometry changes affect the spacing and widths of the wire resonances. In the case of 
odd-numbered wires the transmission is on-resonance, and hardly affected by the resonance shapes, 
whereas for even-numbered wires the transmission is off-resonance and sensitive to the resonance 
shapes. Predicting the amplitude of the conductance oscillation requires a first-principles calculation 
based upon a realistic structure of the wire and the leads. A simple tight-binding model is introduced 
to clarify these results. 
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I. INTRODUCTION 

Recent progress in fabricating conductors of atomic di- 
mensions has stimulated a large number of experimental 
and theoretical studies on these nanoscale devicesi Con- 
ductors whose cross section contains only a small number 
of atoms are commonly called "atomic wires" . Clear evi- 
dence that the fundamental limit of a one atom cross sec- 
tion can be reached, has been presented for gold atomic 
wiresJi^ Over the last decade the electronic transport in 
atomic wires made of various metals has been character- 
ized in great detail experimentally^'S'SiS^ Such wires 
have conductances of the order of the quantum unit 
Go = e^ /irh, so the description of their transport proper- 
ties, as well as of their atomic and electronic structures, 
requires a full quantum-mechanical treatment^ 

Simple theoretical schemes have been proposed, in 
which the atomic wire is described by a jellium^^ or a 
tight-binding model. ^^ At present a first-principles ap- 
proach based on density functional theory (DFT) gives 
the most advanced description of the geometry and elec- 
tronic structure of atomic wires. Several theoretical 
methods that combine the Landauer-Biittiker formal- 
ism with DFT, have been developed in order to solve 
the quantum transport problem in terms of scattering 
amplitudes>i2*i^4i^*iSiiLi§: Alternatively, a Green's func- 
tion formalism is commonly used for solving the trans- 
port problemii"i^"i^'i^^i^'^i24i^'^i2»i^7,2«i2«,^nr-^i Both these 
approaches are in fact completely equivalent in the case 
of nonintcracting electrons.'^ 

Atomic wires that have a cross section of just one atom, 
so-called "monatomic" wires, are the ultimate examples 
of quasi-one-dimensional systems. Here the effects of 
a reduced dimensionality are expected to be most pro- 



nounced. A priori the existence of monatomic wires is 
not obvious. Such free-standing one-dimensional struc- 
tures might be unstable because of the low coordination 
number of the atoms in the wire. Molecular dynamics 
simulations based upon an effective medium model^ or a 
tight-binding modelS^i^ have been used to study the sta- 
bility of a wire as a function of its elongation. Since such 
simulations use highly simplified interatomic potentials, 
they aim at providing a qualitative understanding of the 
wire formation. A more quantitative description can be 
provided by first-principles DFT calculationsr^*^Si22i^2i2S 
but then only relatively small systems can be handled. A 
multitude of different structures has been studied, such 
as dimerized, zigzag, and helical wires 42*^i*^2il2ii4i^ 

One of the most striking features of monatomic wires is 
the nonmonotonic behavior of the conductance as a func- 
tion of the number of atoms in the wire4Siil Such a be- 
havior has been predicted by Lang for wires consisting of 
monovalent atoms4i His model assumes a chain of atoms 
suspended between two planar semi-infinite jellium elec- 
trodes. The conductance predicted by this model is much 
lower than the quantum unit, which disagrees with ex- 
periments on monovalent atomic chains. ^'^'^ However, 
the model can be modified in a simple way by adding 
a basis consisting of three atoms on top of the jellium 
electrodeSf^-^ This reduces the charge transfer between 
the wire and leads and it reduces spurious reflections at 
the wire-electrode interfaces. The conductance of a one- 
atom wire is then close to the quantum unit, in agreement 
with experiment. 

From the Friedel sum rule it can be shown that the 
conductance of an atomic chain consisting of monova- 
lent atoms exhibits a regular oscillation with respect to 
the number of atoms in the chain:— Moreover, assuming 



mirror reflection and time-reversal symmetries together 
with local charge neutrality of the wire, the conductance 
of wires with an odd number of atoms is expected to 
be very close to the quantum unit.'^" The period of the 
conductance oscillation then equals two atoms. This has 
been confirmed in conductance calculations for wires con- 
nected to jellium leads via an atomic basis ^ 

The oscillating behavior of the conductance as a func- 
tion of the wire length has been observed experimen- 
tally in gold wires, ^ and a small parity effect has re- 
cently been found in silver wires. ^^ Oscillations have also 
been observed in wires consisting of atoms with a higher 
valency iSiS 

Since the discovery of the parity effect in wires of 
monovalent atoms there has been a discussion on how 
sensitive the conductance oscillation is to the geome- 
try of the wire and the contacts. The general argu- 
ments given in Ref. |43 suggest that the conductance 
for odd-numbered wires should always be higher than 
for even- numbered ones, provided that the wires are lo- 
cally charge neutral. For a sodium wire connected to 
(artificial) fee sodium electrodes the charge transfer has 
been estimated and it is found to be rather smalli^ Con- 
ductance calculations based upon a scattering approach 
have recently been performed for short sodium wires at- 
tached to sodium electrodes with a more realistic bee 
structure>i2iS The phase of the conductance oscillation 
obtained in these calculations, is consistent with that 
found in the jellium electrode calculations i^2i^ 

Other studies predict however that the conductance 
oscillation found in wires of monovalent atoms is very 
sensitive to the geometry. Even the phase of the oscilla- 
tion can be reversed such, that even-numbered wires have 
a larger conductance than odd-numbered onesi^i*^SiSiS. 
Using the Friedel sum rule to calculate the conductance 
of a wire connected to jellium electrodes it has been found 
that the phase of the conductance oscillation is reversed 
if the jellium leads become sufficiently sharp i^ This has 
not been confirmed by later calculations using a scat- 
tering approach to calculate the conductance, which give 
results that are consistent with Lang's findings for planar 
jellium electrodesiSMl 

From strictly one-dimensional linear combination of 
atomic orbitals DFT (LCAO-DFT) calculations, i.e. 
sodium chains coupled to one-dimensional metallic leads, 
it has been argued that there is a critical distance be- 
tween the wire and the leads where the conductance 
oscillation changes its phase and even-numbered chains 
become more conductive than odd-numbered ones.^^ A 
change of phase has also been predicted to occur upon 
elongating the wire by adding atoms. In short wires the 
odd-numbered chains then have the higher conductance 
and in long wires the even-numbered ones have the higher 
conductance. Similar effects have also been claimed re- 
cently for particular atomic configurations in calculations 
using three-dimensional leadsi^i*^ 

In conclusion, some of the results that appeared in the 
literature regarding the even-odd conductance oscillation 



in monatomic wires seem to be contradictory. In this pa- 
per we present the results of conductance calculations for 
monatomic sodium wires in order to investigate the ef- 
fect of the wire geometry and the wire-lead coupling on 
the conductance oscillation. Since sodium has a simple 
electronic structure, a sodium wire is one of the simplest 
examples of an atomic wire. As such it is an important 
reference system for studying wires with a more compli- 
cated electronic structure, and it can be used as a system 
for benchmarking theoretical and computational tech- 
niques. We perform first-principles conductance calcula- 
tions based on the mode-matching technique^'' on sodium 
wires suspended between sodium electrodes, while sys- 
tematically varying the atomic configuration of the wire 
and that of the wire-lead contacts. The entire system 
consisting of the wire and the semi-infinite electrodes, is 
treated fully atomistically. 

We find that the parity effect, i.e. the even-odd con- 
ductance oscillation, is very robust with respect to chang- 
ing the structure of the wire and to varying the strength 
of the coupling between the wire and the leads. In the 
conductance of long wires we find no tendency to a phase 
change in the even-odd oscillation. The conductance is 
analyzed using the electronic levels of free-standing wires 
in order to interpret the parity effect in terms of trans- 
mission resonances. In addition, we analyze our first- 
principles results using a simple tight-binding model. In 
particular, we show that local charge neutrality of the 
sodium wires provides a strong constraint on the phase 
of the conductance oscillation for all atomic structures 
considered. In absence of a significant charge transfer 
between the wire and the leads, a transmission resonance 
is pinned at the Fermi energy for wires containing an odd 
number of atoms, which leads to a conductance close to 
one quantum unit. Obtaining quantitative values for the 
conductance, particularly for even-numbered wires, re- 
quires well-converged first-principles calculations using a 
realistic structure of the wire and the leads. 

The structure of this paper is as follows. In Sec.Hllwe 
discuss the geometry of infinite and finite sodium wires. 
The even-odd oscillation of the conductance is discussed 
in general terms in Sec. lIIII We investigate the effects on 
the conductance of varying the wire geometry and the 
contacts between wire and leads in Sec. IIVI Current- 
voltage characteristics are analyzed in Sec.|3 In Sec. I VII 
we compare our results to those obtained in previous 
studies. A summary and conclusions are presented in 
the last section. The important technical detail of k- 
point sampling is discussed in Appendix IXI 



II. STRUCTURE OF SODIUM WIRES 

In this section we investigate possible structures of 
sodium monatomic wires by DFT total energy calcula- 
tions in combination with geometry optimizations. DFT 
total energies are calculated with the PW91 general- 
ized gradient approximation (GGA) functional^^ and the 



TABLE I: Na-Na nearest neighbor bond length (in ao) for 
sodium wires, bulk sodium and the sodium dimer, compared 
to all-electron calculations and to experiment. For the zigzag 
chain also the bond angle is given. 
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FIG. 1: (Color online) Supercell representing a two atom 
sodium wire between two sodium leads terminated by (001) 
surfaces. The wire is connected to electrodes via an apex 
atom placed on top of each electrode in a hollow site. 



projector augmented wave (PAW) method, ^°-'-l' as im- 
plemented in the Vienna Ab initio Simulation Package 
(VASP).SiSigi We use a standard frozen core PAW po- 
tential and a plane wave basis set with a kinetic energy 
cutoff of 24 Ry. A Methfessel-Paxton smearing is applied 
in integrations over the Brillioun zone with a smearing 
parameter a = 0.1 eV. First we discuss the structure of 
the infinite wire, and then that of finite wires connected 
to bcc sodium electrodes. 



A. Infinite wires 

An orthorhombic supercell is used with cell parame- 
ters perpendicular to the wire direction equal to 17 A. 
Parallel to the wire the cell parameter is optimized using 
24 k-points to sample the Brillouin zone along the wire 
direction for a cell containing two atoms. 

First we consider a linear wire geometry; the optimized 
Na-Na bond length is given in Table |l| We have checked 
that a symmetry breaking in the form of a Peierls dis- 
tortion, is negligible for Na-Na bond lengths near the 
equilibrium interatomic distance, which agrees with cal- 
culations on other monatomic wires Ji^ Only for sodium 
wires that are stretched to interatomic distances larger 
than 6.91ao a Peierls distortion takes place, accompa- 
nied by a metal-insulator transition. In order to test the 
accuracy of the calculations we have also calculated the 
optimized Na-Na bond length in bulk sodium and in the 
sodium dimer. The accuracy is found to be better than 
1% as compared to the experimental values, see Table Q] 

One-dimensional chains are often unstable with respect 
to a deformation in the transverse direction, which results 
in a zigzag structure. A linear conformation is preferred if 
the interatomic distance exceeds a critical value^i Indeed 
we find that a sodium wire with a zigzag structure has 
a lower energy; its geometry is presented in Table ^ It 
is in reasonable agreement with the one obtained in a 
previous DFT-LDA calculation^^. A bond angle of ~ 60° 
is typically found also in other monatomic wires4^ Upon 
stretching the wire a transformation from a zigzag to a 
linear geometry takes place as soon as the interatomic 
distance in the linear wire becomes > 6. Sap. In Ref. |42 a 
zigzag structure has been found in gold and copper wires, 
whereas in potassium and calcium wires it exists only 



under compression. The stability of the zigzag geometry 
has, therefore, been related to the presence of directional 
d-bonds in gold and copper. However, our results show 
that a sodium wire behaves similarly, which suggests that 
the occurrence of a zigzag geometry is not a result of d- 
bonds only. 

Whereas the lowest energy structure is paramagnetic, 
in Ref. ;6& two additional local minima have been found 
corresponding to ferromagnetic structures with magnetic 
moments ~ 0.02 /ib and 0.12 /iB, respectively. Our lowest 
energy (zigzag and linear) structures are always param- 
agnetic. A magnetic ordering occurs for zigzag structures 
when the Fermi level crosses two energy bands instead of 
just one band, but the magnetism disappears rapidly as 
soon as the wire is stretched sufficiently. Small magnetic 
moments have also been found in calculations on gold 
wires. ^^ No trace of magnetic effects has been observed 
in recent conductance measurements performed in mag- 
netic fieldsi^ Since both experiment and theory favor 
nonmagnetic structures, we will only consider nonmag- 
netic sodium wires in the following. 



B. Finite wires 

In this section we discuss the structure of a finite 
monatomic sodium wire suspended between two elec- 
trodes. A reasonable approach is to study the structure 
of the wire near its equilibrium geometry, which corre- 
sponds to the most stable chemical bonding. We use 
the equilibrium geometry of the infinite linear wire as a 
starting point for finite wires. The electrodes consist of 
bulk Na in the (001) orientation. To calculate the struc- 
ture we use a periodic supercell that consists of a slab of 
five layers of sodium for the electrodes. On top of each 
electrode surface an apex atom is placed in a hollow site 
and a linear wire bridges the two apex atoms as shown 
in Fig. n We use al2xl2x4 k-point sampling of the 
supercell. During the geometry optimization the atoms 
in the wire, the apex atoms, and the atoms in the top 
surface layer are allowed to relax. 

The results of the geometry optimization for wires of 
different lengths are given in Tabled We will discuss the 
most prominent features of the wire geometries starting 
from the electrodes. All structures have mirror symme- 
try with respect to a plane through the center of the wire. 



TABLE II; Optimized bond lengths (in ao) for a Na wire 
suspended between Na electrodes. The columns indicate the 
number of atoms in a wire. The row labels i-j indicate the 
distance between the i'th and j'th atom in a wire; d is the 
average bond length. A-1 indicates the distance between the 
apex atom and the wire, Ll-A the distance between the apex 
atom and the surface layer, and L1-L2 the distance between 
the two top layers of the electrode. The in-plane Na-Na dis- 
tance in the top layer is given in the bottom row. 



bonds 



6 



7 8 



9 



1-2 6.35 6.45 6.34 6.41 6.34 6.41 6.33 6.37 

2-3 6.45 6.46 6.45 6.45 6.43 6.44 6.44 

3-4 6.34 6.45 6.36 6.41 6.37 6.40 

4-5 6.41 6.45 6.41 6.45 6.42 

5-6 6.34 6.43 6.37 6.42 

6-7 6.41 6.44 6.40 

7-8 6.33 6.44 

8-9 6.37 

d 6.35 6.45 6.38 6.43 6.39 6.42 6.39 6.41 

A-1 6.53 6.59 6.51 6.54 6.49 6.51 6.47 6.48 6.46 

Ll-A 3.72 3.64 3.60 3.58 3.52 3.52 3.45 3.43 3.38 

L2-L1 4.03 4.03 4.02 4.03 4.02 4.03 4.02 4.06 4.02 

in-plane 8.17 8.19 8.20 8.23 8.23 8.25 8.26 8.26 8.31 



parallel to the electrode surface. We emphasize that this 
symmetry is not forced upon the system, but is the result 
of the geometry optimization. The top layer of the elec- 
trode relaxes slightly outwards; the distance between the 
top two layers, Ll-L2^ 4.03ao, is somewhat larger than 
the bulk value 3.99ao. The distance between the apex 
atom and the surface Ll-A decreases with the length of 
the wire, which indicates a growing bond strength. The 
distance between the apex atom and the first atom of 
the wire A-1 is always larger than the maximum bond 
length between atoms in the wire. This indicates that 
bonding within the wire is stronger than bonding to the 
electrodes. For gold wires the opposite has been found, 
i.e. the A-1 distance is shorter than the average bond 
length 1^ 

Focusing upon the interatomic distances between 
atoms in the wire. Table E] clearly shows that even- 
numbered wires exhibit dimerization, i.e. an alternation 
between short and long bonds. A similar tendency is 
found in odd-numbered chains, but they have a topo- 
logical defect, i.e. a "kink", in the center of the wire. 
The average bond length of ^ 6.40ao in even and in odd- 
numbered wires is larger than the optimized bond length 
of 6.30ao in the infinite wire. The infinite chain does 
not show a dimerization until the average bond length 
is larger than ^ 6.91ao (see the previous section). This 
strongly suggests that dimerization in finite chains is en- 
forced by their boundaries. A qualitatively similar be- 
havior has also been found in finite gold wiresi^ 

Since dimerization of a finite wire is associated with 
its bonding to the electrodes, one needs to check how 



sensitive the dimerization pattern is to the connection 
between wire and electrode. We have performed calcula- 
tions on larger lateral supercells, and have also made the 
connection more graduate by putting a base of four atoms 
between the electrode surface and the apex atom. These 
structural variations give essentially the same bonding 
pattern in the wires, i.e. even-numbered wires are dimer- 
ized, and odd-numbered wires additionally have a kink 
in the center. The distance A-1 between apex atom and 
wire stays larger than the interatomic distances in the 
wire. These distances can be modified by stretching or 
compressing the wire, but the dimerization pattern is ro- 
bust. In conclusion, the optimized wire structures pre- 
sented in Table ITU can be considered as typical structures 
that are formed by sodium finite chains suspended be- 
tween two semi-infinite electrodes. 



III. CONDUCTANCE OSCILLATION 

Our calculations of the conductance are based on the 
mode-matching technique and we use a real-space finite- 
difference representation of the Kohn-Sham Hamiltonian 
and the wave functions Ji As a first step, the one-electron 
self-consistent potentials of the bulk leads and the scat- 
tering region containing the wire are obtained from DFT 
calculations. Subsequently the scattering problem is 
solved at the Fermi energy by matching the modes in 
the leads to the wave function in the scattering region. 
The conductance G can be expressed in terms of normal- 
ized transmission amplitudes i„_„' using the Landauer- 
Biittiker iorumlaS^ 



G = Go V \t 



/ ^ I'^n,?!' 



(1) 



where n and r\l label the right-going modes in the left and 
right leads, respectively and Gq = e^ jirh. An efficient 
implementation of a high-order finite-difference scheme 
for solving the scattering problem is discussed in Ref . [l^- 

In more detail, the one-electron potentials of the leads 
and the scattering regions are extracted from two DFT 
calculations for bulk bcc sodium and for the supercell 
shown in Fig. [3 respectively. For these calculations we 
use a local TrouUier-Martins pseudopotential2^ with a 
core radius r^ — 2.95 oq; only the 3s electrons of sodium 
are treated as valence electrons. All plane waves are in- 
cluded up to a kinetic energy cutoff of 16 Ry. We use 32"^ 
k-points to sample the Brillouin zone (BZ) of the cubic 
hcc unit cell of bulk sodium. In our supercell calcula- 
tions, S'^ k-points are used to sample the lateral BZ in 
case of a 2 X 2 supercell, and 6^ k-points in case of a 3 x 3 
supercell. In all cases a Methfessel-Paxton smearing with 
a — 0.1 eV is applied. Total energies are converged to 
within 5 x 10~^ Hartrees. 

One assumes that the leads outside the scattering re- 
gion are perfectly crystalline bulk material. So at the 
edges of the scattering region, the potential should join 




FIG. 2: (Color online) Structure of an atomic wire consisting 
of two sodium atoms between two sodium leads terminated by 
(001) surfaces. The boundaries of the supercell are indicated 
by dashed lines. Bulk atoms are indicated by yellow (light 
grey) balls and atoms in the scattering region by blue (dark 
grey) balls, respectively. 
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FIG. 3: Conductance (in units of Go) as a function of the 
number of atoms in the atomic chain. All atomic bond 
lengthes in the system are equal to the bulk value d = 6.91ao. 



smoothly to the potentials of the bulk leads. We have 
checked that this is the case. Enlarging the scattering 
region by including two extra atomic layers in each lead 
changes the results reported for the conductance only by 
~ 1.5% for even-numbered wires and < 0.5% for odd- 
numbered wires. The Fermi energy is extracted from 
the bulk calculation.^^ The only parameters in calcu- 
lating the conductance within the mode-matching finite- 
difference scheme are the order N of the finite-difference 
approximation of the kinetic energy (i.e., the second 
derivative) and the spacing h^^y^z between the real-space 
grid points. We use A'' = 4 and hx,y,z = O.SOoo; for de- 
tails and convergence tests we refer to Ref. [T^- The total 
transmission is averaged over the k|| -point grid of the lat- 
eral BZ of the supercell. To calculate the transmission it 
is important to apply a proper k|| -point sampling. This 
will be discussed in Appendix fXl Most calculations are 
done for a 2 X 2 lateral supercell. Enlarging the super- 
cell changes the conductance only marginally as will be 
discussed in Sec. IIVBI 



The electron transport in the crystalline leads is bal- 
listic, i.e. an electron goes through the leads without 
any scattering. The transport properties of a monatomic 
wire suspended between two leads depend upon three 
factors; the number of atoms in the chain, the geometry 
of the wire, and the contact between wire and leads. In 
Sec. IIVI we will discuss how these factors influence the 
conductance. In the present section we will analyze the 
conductance of monatomic sodium wires in a reference 
geometry, where all Na-Na bond lengths are chosen to 
be equal to the bulk value 6.91ao, see Table 13 As in 
the previous section we attach a finite atomic wire to the 
leads via two apex atoms, which then have a coordina- 
tion number 5. All atoms in the wire have a coordination 
number 2. 



A. First-principles calculations 

The calculated conductance as a function of the num- 
ber of atoms in the atomic chain is given in Fig. |3| Since 
a sodium atom has valence one, both the infinite sodium 
chain and bulk sodium have a half-filled band, and the 
infinite wire has one conducting channel at the Fermi 
leveli^ The conductance of the infinite chain is then 
equal to the quantum unit Go, and the conductance of 
finite wires is < Go- As can be observed in Fig. Othe 
conductance exhibits a regular oscillation as a function 
of the number of atoms in the wire. The conductance is 
very close to Go for odd-numbered wires, and for even- 
numbered wires it is '^ 10% lower. Such a behavior of 
the conductance in atomic-sized conductors is very dif- 
ferent from ohmic behavior in macroscopic conductors; it 
expresses the quantum nature of the electron transport 
at the nanoscale. 

In order to interpret the even-odd oscillation we have 
calculated the conductance as a function of energy for 
wires of different length. The results for monatomic wires 
consisting of four and five atoms are shown in Fig.^ Res- 
onant peaks in the conductance can be clearly identified. 
Qualitatively they correspond to energy levels of a free- 
standing Na wire, which are shifted and broadened into 
resonances by the interaction of the wire with the leads. 
To illustrate this, the calculated energy levels of free- 
standing wires of four and five atoms are shown as bars 
in Fig. ^ The levels are sufficiently close to the resonant 
energies to warrant an interpretation of the conductance 
in terms of a transmission through levels of the wire. As 
is clearly observed in Fig. 01 the Fermi level is in between 
two resonant peaks for a four atom wire and right on 
top of a resonance for a five atom wire. By calculating 
the conductance as a function of energy for wires of dif- 
ferent length it can be shown that this observation can 
be generalized. The Fermi level is between resonances 
for even-numbered wires and on top of a resonance for 
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FIG. 5: Energy levels for odd- and even- numbered wires. The 
Fermi level is in the middle of the HOMO-LUMO gap for 
even-numbered wires, and it aligns with the HOMO level for 
odd-numbered wires. 
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FIG. 6: Tight-binding representation of the n-atomic wire 
attached to two semi-infinite one- dimensional leads. 



FIG. 4: (Color online) Conductance (in units of Go) as a 
function of energy for monatomic wires consisting of four (top 
figure) and five (bottom figure) atoms. The red (grey) bars 
correspond to the energy levels of free-standing wires. _E = 
corresponds to the Fermi level. 



odd-numbered wires. 

An intuitive picture of the transmission through the 
energy levels is then presented by Fig. |31 Odd-numbered 
wires have a highest occupied molecular orbital (HOMO) 
that is half filled. Perfect transmission through this state 
takes place if the Fermi level aligns with the HOMO. In 
even-numbered wires the HOMO is completely filled and 
separated by a gap from the LUMO (lowest unoccupied 
molecular orbital) level. The Fermi level is then in the 
HOMO-LUMO gap. The position of the Fermi level with 
respect to the HOMO level causes the oflf and on resonant 
behavior of the conductance as a function of the wire 
length, which is causing a regular even-odd oscillation of 
the conductance. In the next section we will study this 
intuitive model in somewhat more detail by means of a 
simple tight-binding model. 



B. Tight-binding model 

To support the intuitive picture presented in the previ- 
ous section we use a simple tight-binding model as shown 
in Fig. in which the leads are modeled as quasi-one- 
dimensional systems described by effective parameters. 



Here eoi P are the on-site energies and nearest neighbor 
hopping coefficients of the leads, and eg, /3' are the corre- 
sponding parameters of the wire. The coupling between 
the left (right) electrode and the atomic chain is given by 
the hopping coefficient /3c (/3c)- 

If the system has mirror symmetry, the coupling is 
symmetric, i.e. fie = /3c- The leads and the chain are 
made of the same material (sodium). If one assumes that 
all atoms are neutral (local charge neutrality), then it is 
not unreasonable to set Eq = Eq- The conductance can 
be calculated analytically for this model by the mode- 
matching technique. '^^ The modes can be labeled by a 
wave number fc in ID Brillouin zone of the leads. The fa- 
miliar relation £^ = 60 + 2/3cos(A;a) gives for a half-filled 
band the Fermi energy Ep = Eq and the Fermi wave 
number kp = 7r/2a. The parameter /3 can be used as a 
scaling parameter. In the following all energy parameters 
eQ,e'Q,P', f3c, /3'c are in units of (3. The conductance of a 
wire at the Fermi energy consisting of n atoms is given 
by 
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The conductance for odd-numbered wires is equal to the 
quantum unit, and it is smaller than the quantum unit for 
even- numbered wires (unless /3^ = f3'). This corresponds 
to the situation shown in Fig. 

It is instructive to study some other consequences of 



the tight-binding model. If Ae = Eq — e'q 7^ then a 
charge transfer will take place between the leads and the 
wire. The conductance calculated at the Fermi energy 
for a one-site wire (n = 1) and a two-site wire [n = 2) 
become, respectively, 
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According to Eq. (j3Jl a nonzero Ae suppresses the trans- 
mission through a one-site wire. The transmission is 
shifted "off resonance" and the conductance becomes 
smaller than the quantum unit. However, the coupling 
between wire and lead also causes a broadening of the res- 
onance, which is proportional to [3c- This broadening par- 
tially compensates for the decrease of the conductance. 
If the coupling is sufficiently strong, i.e. 4/3^ ^ Ae^, 
then the conductance is again close to the quantum unit. 
In the limit of weak coupling, i.e. 4/3^ ^ Ae^, the con- 
ductance goes to zero with decreasing (3c for any nonzero 
Ae. The conductance as a function of Pc is shown in 
Figs.Ua) and[7|[b) for two different values of Ae. 

The conductance of a two-site wire, see Eq. (0}, be- 
haves qualitatively different as a function of the coupling 
strength (3c- In the weak coupling limit, i.e. (3c ^ 
{(3' ± Ae)^, corresponding to /3c < 1 in Fig. CJa), the 
conductance goes to zero with decreasing (3c and the de- 
crease is faster than for a one-site wire. Note that this 
only holds for Ae <C /3'. If Ae ^ (3' then the conduc- 
tance decreases more slowly with decreasing /3c for a two- 
site wire than for a one-site wire, see the range (3c < I 
in Fig. Efb). If the couphng between wire and lead is 
strong, i.e. /?c ^ (/3' ± Ae)^, corresponding to /3c > 1 
in Figs. [3a) and|7{b), then the conductance always de- 
creases with increasing /3c. This is due to a phenomenon 
called "pair annihilation" of resonances,™ which happens 
if the resonance widths become larger than the spacing 
between the resonances. In a one-site wire this cannot 
happen, since there is only one resonance. Between the 
strong and weak coupling regimes there is a value of (3c 
(close to 1) where the conductance of a two-site wire is 
equal to the quantum unit, see Figs.[7{a) and[7Ib). 

The conductance of longer wires, i.e. n > 2, can be 
interpreted along the same lines. For small Ae, the odd- 
numbered wires resemble the one-site wire and the even- 
numbered wires resemble the two-site wire, as shown in 
Fig. [7Ja) . For a very large range of coupling strengths (3c 
one obtains an even-odd oscillation in the conductance of 
a nearly constant amplitude. The conductance of odd- 
numbered chains is close to the quantum unit and that 
of even-numbered chains is smaller by an amount that 
depends upon the coupling between wire and lead. Ap- 
parently, this is the case that corresponds to the results 
of our first-principles calculations, see Fig. |31 

If Ae becomes larger, the conductance of all wires as 
a function of /3c becomes qualitatively similar to that of 




FIG. 7: (a) Conductance (in units of Go) at the Fermi energy 
as a function of /?c (in units of /?) for wires with n = 1 (thick 
solid line), 2 (large dots), 3 (thin solid), 4 (thin dashed), and 
5 (thick dashed) atoms; Ae = 0.02; /3' = 1, (b) As (a) with 
Ae = 0.7; /?' = 1. 



the two-site wire, see Fig. |7{b) (except the one-site wire, 
of course). The amplitude and even the phase of the 
conductance oscillation as a function of the wire length 
then strongly depends upon the coupling /3c of the wire 
to the lead. For instance, if /3c < 0.7 in Fig. [^b), the 
conductance of even-numbered wires is higher than that 
of odd-numbered wires and all conductances are smaller 
than the quantum unit. Note that if Ae is significant, it 
will be accompanied by a significant charge transfer be- 
tween wire and leads. Whether this situation occurs can 
be studied by self-consistent first-principles calculations. 



IV. STABILITY OF CONDUCTANCE 
OSCILLATION 

First-principles calculations on a reference geometry 
give a regular even-odd oscillation of the conductance as 
a function of the wire length, as discussed in the previ- 
ous section. The odd-numbered wires have the highest 
conductance, close to the quantum unit Gq. The sim- 
plified tight-binding model suggests that the geometry 
might influence the amplitude and even the phase of the 
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FIG. 8: Conductance (in units of Go) as a function of the 
number of atoms in the wire. The triangles correspond to 
stretched wires with d — 6.91ao, squares to compressed wires 
with d — 5.82ao, and diamonds to wires with d — 6.30ao. 
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FIG. 9: Conductance (in units of Go) as a function of the 
number of atoms in the wire. The triangles correspond to 
wires with d — 5.82ao and dc — 6.91ao and the squares to 
wires with d = 6.91ao and dc = 5.82ao. 



conductance oscillation. In Sec. E] we have shown that 
monatomic sodium wires can have a linear, zigzag, or 
dimerizcd geometry, depending upon the boundary con- 
ditions. In this section we study the influence of the wire 
geometry upon the conductance. In particular, we focus 
on the question of whether the phase of the even-odd 
conductance oscillation is robust to modifications of the 
wire geometry. 



A. Tension or compression of linear wires 

We consider a uniform tension or compression of the 
wire. As in the previous section the Na-Na distance be- 
tween atoms in the wire is kept at a uniform value d and 
the distance dc between the apex atom and the wire is 
equal to d. As reference we use the results shown in Fig.O 
where d — 6.91ao, which corresponds to the bond length 
in bulk sodium. If we take the equilibrium bond length 
d — 6.30ao of the infinite linear wire as a characteristic 
bond length, then d = 6.91ao corresponds to a wire under 
tension, i.e. a stretched wire. A wire with d — 5.82ao, 
which corresponds to the equilibrium bond length of a 
Na2 molecule, is a wire under compression. 

The results of first principles calculations of the con- 
ductance as a function of the wire length for d — 
5.82ao,6.30ao,6.91ao are presented in Fig.|S| In all cases 
the conductance exhibits a regular even-odd oscillation 
and the conductance of the odd-numbered wires is close 
to the quantum unit. The conduction of the even- 
numbered wires is smaller than the quantum unit and it 
depends only weakly on the number of atoms in the wire. 
According to the tight-binding analysis in Scc. lIIIBl this 
suggests that the charge transfer between wire and lead, 
represented by Ae in Eqs. ^ and Q), is very small, see 
Fig. Ha). 

The amplitude of the oscillation decreases with de- 
creasing d. Two opposing effects influence the conduc- 



tance if we decrease the interatomic spacing d in the 
wire. Firstly, the spacing between the resonant levels 
increases. In tight-binding terms the parameter /?' in- 
creases, which tends to decrease the conductance of even- 
numbered wires, see Eq. 0. Secondly, since dc = d, the 
resonances become broader if we decrease the distance 
between the wire and the lead. Again in tight-binding 
terms the parameter f3c increases, which tends to increase 
the conductance of even-numbered wires, cf. Eq. Q. — 
According to Fig. [S] the effect of the resonance broad- 
ening upon the conductance is larger than the effect of 
increased resonance spacing. 

One can strengthen this analysis by varying the in- 
teratomic distance d in the wire, and the distance dc 
between the apex atom and the wire independently. Fig- 
ure El shows the calculated conductance for a wire with 
d < dc, i.e. d = 5.82ao and dc = 6.91ao. The conduc- 
tance oscillations are quite large, which can be attributed 
to the increased resonance spacing discussed in the pre- 
vious paragraph. A small d results in a large spacing 
between the resonant levels of the wire. Therefore, the 
transmission of even-numbered wires, which is off reso- 
nance at the Fermi level, is low, whereas the transmission 
of odd-numbered wires stays on resonance and is high. 

If we calculate the conductance for a wire with d > dc, 
i.e. d = 6.91ao and dc — 5.82ao, we see in Fig.lHlthat the 
conductance oscillation is strongly suppressed. It can be 
attributed to the resonance broadening. If the coupling 
between the wire and the lead is strong, the resonances of 
the wire are wide. The transmission in even-numbered 
wires is then relatively high, whereas the transmission 
in odd- numbered wires stays close to the quantum unit. 
Figure El shows that in the case of a strong coupling be- 
tween wire and lead the amplitude of the even-odd os- 
cillation in the conductance can become very small. Ac- 
cording to the tight-binding model, Eq. ||21l, this hap- 
pens if Pc/0''^ ~ 1- Note that such a strong couphng 
is less likely for sodium monatomic wires with optimized 
geometries, because the results discussed in Sec. HTIindi- 
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FIG. 10: (Color online) Structure of an atomic wire consisting 
of two sodium atoms between two sodium leads terminated by 
(001) surfaces. The atomic wire is connected to each surface 
via a five atom pyramid. The boundaries of the 3x3 supercell 
are indicated by dashed lines. Bulk atoms are indicated by 
yellow (light grey) balls and atoms in the scattering region by 
blue (dark grey) balls, respectively. 
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FIG. 11: Conductance (in units of Go) as a function of the 
number of atoms in the atomic wire. Squares refer to the 3x3 
supercell with five atom pyramid contacts, diamonds to the 
3x3 supercell with one apex atom contacts, and triangles to 
the 2x2 supercell with one apex atom contacts. 



cate that d < dc- According to Ref. |33 d > dc in gold 
monatomic chains, which might explain the small ampli- 
tude of the conductance oscillation found experimentally 
in gold wires 1^. 

In conclusion, stretching or compressing the wire 
changes the amplitude of the conductance oscillation, but 
it preserves its phase and the value of the conductance 
for odd- numbered wires, which is close to unity. 



B. Contact geometry 

The coupling between the wire and the leads could be 
influenced by the detailed geometry of the two contacts 
between the wire and the leads. *^ Since the geometries of 
the wire-lead contacts are not known from experiment, 
it makes sense to study the sensitivity of the calculated 
conductances to these geometries. So far in our calcula- 
tions we have modeled both contacts by one apex atom 
placed in a hollow site on the (001) surface in a 2 x 2 
lateral supercell. A more "gradual" contact is formed by 



FIG. 12: (a) Conductance (in units of Go) as a function of the 
number of atoms in the atomic wire with interatomic distance 
d — dc = 6.91ao. The triangles (dashed line) correspond to 
"direct" contacts between wire and leads. The squares (solid 
line) correspond to a direct contact at one end of the wire and 
a one apex atom contact at the other end. (b) As (a) with 
d = dc — 5.82ao. 



a five atom pyramid placed on (001) surface, as is shown 
in Fig.^1 This requires using (at least) a 3 x 3 supercell. 
To check that the size of the supercell does not influence 
the results, we have also done calculations for a 3 x 3 
supercell with one apex atom contacts. The calculated 
conductances are shown in Fig. ^] The conductance of 
monatomic sodium wires seems to be relatively insensi- 
tive to the contact geometry. At the same time it shows 
that the results obtained with the 2x2 supercell are 
converged. 

Another way of modifying the contacts is to remove 
the apex atoms and position the first and the last atom 
of the wire on top of an atom in the (001) surface layer. 
The calculated conductances are shown in Figs. I12r a') 
and ll2r b') for the interatomic distances d ^ dc = 6.91ao 
and d = dc — 5.82ao, respectively. The results obtained 
with these "direct" wire-surface contacts look very sim- 
ilar to the ones obtained with one apex atom contacts, 
see Fig. |S1 The amplitude of the even-odd oscillation 
is somewhat smaller for the "direct" coupling. Accord- 
ing to the analysis presented in Sec. lIVAl this indicates 
a stronger coupling between wire and leads, or in tight- 
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FIG. 13: Conductance (in units of Go) as a function of the 
number of atoms in the atomic wire. Dashed hue and filled 
rectangles correspond to zigzag wires; straight line and filled 
squares to linear wires with d — 6.91ao. 



FIG. 14: Conductance (in units of Go) as a function of the 
number of atoms in the atomic wire. The squares correspond 
to dimerized wires with alternating bond lengths of 6.91ao 
and 5.82ao; as a reference, the triangles correspond to wires 
with a uniform bond length 6.91ao. 



binding terms, a larger f3c, of. Eq. ^. Note that the 
conductance of one-atom wires in Figs. I12r a) and ll2f b') 
is higher than Gq due to direct tunncUng between the 
electrodes. 

One can also break the symmetry and use a direct 
contact between the wire and one of the leads, and a one 
apex atom contact between the wire and the other lead. 
The calculated conductances are given in Figs. ll2f a) and 
I12r b). A comparison with symmetric direct contacts and 
symmetric one apex atom contacts, see Fig.|S| shows that 
the phase of the even-odd conductance oscillation is the 
same and the amplitude is in between that of the two 
symmetric cases. It means that, besides the already men- 
tioned stronger coupling between wire and lead for the 
"direct" contact, this symmetry breaking has little effect 
on the conductance. 

We conclude that varying the geometries of the con- 
tacts between wire and leads does not have a large effect 
on the regular even-odd oscillation of the conductance. 



C. Wire geometry: zigzag wires and dimerization 

In Sec. m we studied the geometry of infinite zigzag 
chains. In principle, a structural zigzag deformation 
could modify the conductance of a finite monatomic 
^lYe^dM,. Figure El shows the calculated conductance 
of a wire with one apex atom contacts and bond lengths 
d = dc = 6. 9 lap to which a zigzag distortion pattern is 
applied with an amplitude corresponding to 15% of the 
bond length. Compared to straight wires, the conduc- 
tance of zigzag even-numbered wires changes by < 3%, 
whereas the conductance of odd- numbered wires is hardly 
affected at all. Such small effects are in line with results 
reported previously42ii& 

According to the results obtained in Sec. |n] finite 
straight wires with equidistant atoms can spontaneously 
break their symmetry by dimerization. The conductance 
of optimized broken symmetry structures is discussed in 



the next section. Here we study the influence of an exces- 
sive symmetry breaking. We apply a regular dimerization 
pattern to the wire, which consists of an alternation be- 
tween long and short bonds with bond lengths d = 6.91ao 
and d = 5.82ao, respectively. Continuing this pattern 
into the contacts this means that even-numbered wires 
have short dc = 5.82ao contacts to both leads, whereas 
odd-numbered wires have one short dc = 5.82ao contact 
and one long dc = 6.91ao contact. The results are shown 
in Fig. m 

This curve can be analyzed using the tight-binding 
model introduced in Sec. IIII Bl Assuming charge neu- 
trality, i.e. Ae = 0, we need to generalize Eq. (0) to the 
case where the coupling /3c to the left and /3^ to the right 
leads are different. In addition, dimerization in the wire 
leads to an alternation of two hopping coefficients /3' and 
/3". The conductance of an n-atomic wire is then given 
by 



G = Gq — — -^ — 7T, nodd 




n even. (5) 



One notices from Eq. © that even the shortest odd- 
numbered wire, n = 1 has a conductance smaller than 
the quantum unit if /3c ^ /3c- This is observed in our first 
principles results, where the conductance of the n = 1 
wire is 0.95Go, see Fig. E] Furthermore, Eq. lO shows 
that the conductance of both even-numbered and odd- 
numbered wires decreases as a function of increasing n 
if /3' ^ (3". Also this is clearly observed in our first 
principles calculations. Dimerization of an infinite wire 
creates a gap in the density of states, so, in general, one 
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expects that the conductance drops as a function of the 
wire length. A decreasing conductance for longer wires 
has been observed experientially for platinum, but its 
nature has not been clarified yetuS 

Finally, although the conductance for even- and odd- 
numbered wires decreases as a function of wire length, 
its even-odd oscillation is preserved. Assuming /3"//3' = 
Pel (i'c = ^ foi' odd-numbered wires and P" / (3' = Pd fi = 
Pel P' = y for even-numbered wires the tight-binding 
model, Eq. 0, can be fitted to the first-principles re- 
sults. This yields the parameter ratio's x « 0.81 and 
y « 0.76 for the odd- and even-numbered dimerized wires 
presented in Fig. [HI 



D. Optimized geometry 




FIG. 15: Conductance (in units of Go) as a function of the 
number of atoms in the atomic wire. The triangles correspond 
to hnear wires with d — 6.30ao and the squares to wires with 
optimized geometry. 



In previous sections we studied the influence of the 
structure of a monatomic wire upon its conductance by 
varying interatomic distances corresponding to values 
ranging from the Na dimer 5.82ao to the Na bulk 6.91ao 
values. In this section we discuss the conductance for 
wires with optimized geometries, which were obtained 
in Sec. |nj Figure El shows the results from the first- 
principles calculations. As a reference, it also shows the 
results for wires with equidistant atoms corresponding to 
the geometry of an infinite wire. 

It can be observed that the amplitude of the conduc- 
tance oscillation for wires with optimized geometry is 
larger than for the reference wires. This results from 
a slight decoupling of the wire from electrodes, since for 
the optimized wires dc > d, see Table ^ and the dis- 
cussion in Sec. II V Bl As can be seen from Fig. ^1 the 
phase of the even-odd conductance oscillation is also ob- 
served for optimized geometries. Moreover, despite the 
presence of a topological defect in the center of the odd- 
numbered wires, the conductance of odd-numbered wires 
is close to the quantum unit. This is actually predicted 
by the tight-binding model for any odd-numbered wire 
whose geometry has mirror symmetry with respect to a 
plane through the center and perpendicular to the wire, 
provided Ae = 0. 

One might expect the conductance for even-numbered 
wires to decrease with the length of the wire, due to the 
effect of dimerization as discussed in the previous section. 
This does not show in Fig. 1151 because the dimerization in 
the optimized geometry is much weaker. Therefore, the 
effect will show up only in wires that are much longer. 



V. BEYOND LINEAR RESPONSE 

In the linear response regime the current and therefore 
the conductance are fully determined by the electrons at 
the Fermi energy. If a finite bias V is applied then the 



current is given by 



I-Ef+V/2 

I^Go T{E, V) dE, 

Jef-V/2 



(6) 



where the transmission coefficient T (_E, V) depends on 
the energy of the electron E and the voltage V. The 
differential conductance is defined as G{V) = dl IdV . 

In this section we discuss some of the consequences of 
a finite bias. We make the approximation T {E, V) « 
T (E), which is valid for a relatively small voltage in 
the limit that the electronic structure of the wire is not 
changed by the voltage. Examples of transmissions as a 
function of energy are given in Fig. ^ for atomic wires 
with the "reference" geometry d — dc — 6.91ao. The 
corresponding I-V curves, calculated from Eq. lO, are 
given in Fig. E| From the /- V curves we calculate the 
differential conductance and the second derivative of the 
conductance, which are also presented in Fig. 1161 

The conductance varies by less than 5% for biases up 
to ~ ±0.2V, which one might call the linear response 
regime. The conductance decreases monotonically for 
odd-numbered wires and it increases monotonically for 
even-numbered wires for biases up to ~ 0.5V. The oscil- 
lating behavior of the conductance at higher biases re- 
sults from the resonant peaks in the transmission. At 
biases larger than ~ 0.5V the non-self-consistent proce- 
dure probably becomes increasingly inaccurate. ^^ 

An important characteristic of the conductance curve 
that can be measured experimentally, is its second deriva- 
tive at the Fermi energy. If the transmission near the 
Fermi energy E-p = is approximated by a polynomial 
function T{E) = T{0) + T{0)" E'^, then the differen- 
tial conductance is G{V) = Go [T(0) + T{0)"V'^li]. In 
Sec, mil we interpreted the even-odd conductance oscilla- 
tion of the conductance in terms of switching between off- 
and on-resonance behavior. If this is true then the sec- 
ond derivative of the conductance d^Gld'^V = r(0)"/2 
must be positive for even-numbered wires and negative 



12 


















1 


- 




1 


" 








A 




~ 






-1 




• 










: 


^ 












\ 
\ 


: 


o 














\ 
\ 


: 


■^ 


-2 
-3 










1 1 


• 

1 


- 



FIG. 16: Current-voltage characteristics for two (thick dashed 
line), three (thick line), four (thin dashed line) and five (thin 
line) atomic wires with d = dc — 6.91ao. The top, middle 
and bottom figures contain I-V curves, the differential con- 
ductance G(V) = dl /dV as a function of V and the second 
derivative of the conductance d^G{0)/d^V as a function of the 
number of atoms, respectively. 



for odd-numbered wires at the Fermi energy, as shown in 
Fig. ^1 At the same time the first derivative is zero. 



VI. DISCUSSION 

In this section we compare our results to those obtained 
in previous studies. The oscillating behavior of the con- 
ductance of monatomic sodium wires was first suggested 
from calculations using planar jellium electrodes 4i The 
conductance of both even- and odd-numbered wires is 
then significantly lower than the quantum unit and the 



conductance of even-numbered wires is larger than that 
of odd-numbered ones. Both of these features are likely 
to be artifacts of using jellium electrodes, since adding 
atomic bases between wire and jellium leads reverses the 
phase of the conductance oscillation and makes the con- 
ductance of odd-numbered wires approach the quantum 
unit.'*^'^^ Calculations using tip-shaped jellium electrodes 
predict that the phase of the conductance oscillation crit- 
ically depends upon the sharpness of the tips^^ although 
this effect is disputed in recent calculationsiSS 

Using atomistic electrodes with a bulk structure our 
calculations show that the conductance has a regular 
even-odd oscillation, in which the conductance of the 
odd-numbered wires in optimized structures is close to 
the quantum unit and that of even-numbered wires is 
approximately 10% lower. Only the latter is modified 
substantially if the geometry of the wire or the contacts 
between wire and electrodes are changed within reason- 
able bounds. Apparently jellium electrodes cause reflec- 
tions of electrons trying to enter the wire, which results in 
the artifacts discussed above. The amplitude of the con- 
ductance oscillation we find, is smaller by ^ 67% than 
that obtained using jellium electrodes plus atomic bases. 
This would indicate that atomic bases do not completely 
remove the reflections caused by the jellium electrodes. 

Calculations based upon one-dimensional metal elec- 
trodes predict that it is possible to change the phase of 
the conductance oscillation by varying the coupling be- 
tween the wire and the electrodes^ We did not observe 
such an effect for three-dimensional atomistic electrodes. 
In a recent calculation it is found that the conductance 
of odd-numbered wires decreases sharply with increas- 
ing wire length, and even the phase of the conductance 
oscillation can be reversed in long wiresi^ This is not 
confirmed by our calculations, where the conductance of 
odd-numbered wires stays close to the quantum unit and 
the phase of the oscillation is stable. Experiments on 
the even-odd conductance oscillation in monovalent gold 
wires do not reveal a decrease of the conductance in odd- 
numbered wires, and the conductance stays close to the 
quantum uniti^ 

Other studies seem to indicate that the phase of the 
even-odd conductance oscillation does not depend very 
sensitively upon the structure of the electrodes, since the 
same phase is observed in calculations using bcc elec- 
trodes oriented in the (111) direction, ^^ and in calcu- 
lations using electrodes with an artificial fee structure.^ 
The amplitude of the oscillation is much more sensitive, 
however. A previous study on sodium wires suspended 
between sodium electrodes gives an amplitude of only 
1%, which is an order of magnitude smaller than what 
we find using similar geometries. ^^ Although we do not 
know what the cause of this difference is, we observe that 
the amplitude of the even-odd conductance oscillation is 
sensitive to the one-electron potential used in solving the 
scattering problem. This potential is obtained from a 
self-consistent electronic structure calculation. Such cal- 
culations frequently use a convergence criterion applied 
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FIG. 17: Conductance (in units of Go) as a function of the 
number of atoms in the atomic wire. Triangles and squares 
correspond to using a potential obtained from a total energy 
calculation with energy convergence criterion set to AE — 
2 X 10~^ and 5 x 10~^ Hartrees, respectively. 



to the total energy. However, since a variational prin- 
ciple does not apply to the one-electron potential, the 
convergence criterion should be much stricter in order to 
converge the potential. 

As an illustration, Fig. El shows the conductance of 
sodium wires calculated from a potential obtained with 
the usual energy convergence criterion, compared to one 
obtained with a stricter energy convergence criterion. We 
have checked that the result does not change anymore 
if the convergence criterion is even made stricter. This 
figure clearly shows that changes in the potential that are 
caused by small charge transfers can markedly influence 
the amplitude of the even-odd conductance oscillation. 



VII. SUMMARY AND CONCLUSIONS 

We have performed first-principles calculations to 
study the stability of even-odd conductance oscillations 
in a sodium monatomic wire with respect to structural 
variations. An infinite sodium wire can be linear and con- 
sists of equidistant atoms or dimers of atoms, or it can 
have a zigzag structure, depending upon the tension or 
compression applied to the wire. The geometry of finite 
sodium wires, suspended between two sodium electrodes, 
is influenced by boundary effects. Wires comprising an 
even or odd number of atoms are dimerized, but odd- 
numbered wires have a topological defect in the center. 

In the linear response regime the conductance is de- 
termined by the electrons at the Fermi energy. The 
conductance of sodium wires shows a distinct even-odd 
oscillation. The odd-numbered wires have a conduc- 
tance close to the quantum unit Gq — e'^/nh and even- 
numbered wires have a lower conductance. This oscilla- 
tion is remarkably robust, as we show by systematically 
varying the structure of the wires and the geometry of 
the contacts between the wires and the electrodes. The 
phase of the oscillation is not affected by these struc- 



tural variations, i.e. odd-numbered wires have a higher 
conductance than even-numbered ones. Moreover, odd- 
numbered wires have a conductance close to the quan- 
tum unit, unless the structural deformation of the wire 
becomes very large and the contact to the left lead is 
markedly different from that to the right lead. The con- 
ductance of even-numbered wires is much more sensitive 
to the wire geometry. Increasing the interatomic dis- 
tances in the wire and/or strengthening the contacts be- 
tween wire and leads increases the conductance of even- 
numbered wires; increasing the asymmetry between the 
interatomic distances or between left and right contacts 
decreases the conductance. 

These results can be interpreted on the basis of res- 
onant transmission. For odd-numbered wires the Fermi 
energy coincides with a resonance in the transmission, 
whereas for even-numbered wires the Fermi energy is be- 
tween two resonances. Changing the geometry of the wire 
or the contacts affects the spacing between the resonances 
and their widths and therefore it affects the conductance 
of even-numbered wires; decreasing the spacing and/or 
increasing the widths increases the conductancei^ Since 
for odd-numbered wires the Fermi level is pinned at a res- 
onance, their conductance is affected much less by chang- 
ing the wire geometry. 

We have formulated a simple tight-binding model to 
analyze these results. It shows that the even-odd con- 
ductance oscillation is stable with respect to structural 
variations, unless the on-site energies for atoms in the 
wire are substantially different from the on-site energies 
of atoms in the leads. Note that a large difference in on- 
site energies is necessarily accompanied by a significant 
charge transfer between the wire and the leads. The re- 
sults of the first principles calculations demonstrate that 
this is not the case. For wires with equidistant atoms that 
have mirror symmetry with respect to a plane perpendic- 
ular to the wire, Eq. Q shows that if all on-site energies 
are identical, the conductance of odd-numbered wires is 
one quantum unit, whereas that of even- numbered wires 
is determined by the ratio of the wire-lead coupling and 
the atom-atom coupling within the wire. Breaking the 
mirror symmetry, Eq. |SJ| shows that the conductance of 
odd- numbered wires becomes smaller than one unit. The 
symmetry breaking has to be large, however, in order to 
have a sizable effect on the conductance. 

We have also calculated the current-voltage character- 
istics of sodium wires in the low bias regime. The differ- 
ential conductance clearly shows a nonmonotonic behav- 
ior. In particular, the second derivative of the conduc- 
tance has an alternating sign as a function of the number 
of atoms in the wire; even-numbered wires have a posi- 
tive second derivative and odd-numbered wires a negative 
one. This effect can be ascribed to the resonant nature 
of the transmission. It could be used to establish the res- 
onant behavior of the even-odd conductance oscillation 
experimentally. 

Comparison to other work shows that simple jellium 
electrodes do not reproduce the even-odd conductance 
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oscillation correctly. Using atomic bases yields the cor- 
rect phase of the oscillation. The conductance of odd- 
numbered wires is rather stable with respect to varying 
the atomic structure, but that of even-numbered wires 
is sensitive to structural details and the quality of the 
one-electron potential. 
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APPENDIX A: k-POINT SAMPLING 

In this Appendix we discuss the effect of ky-point sam- 
pling on the conductance. In modeling a conductor be- 
tween two semi-infinite electrodes one usually assumes 
a supercell geometry in the lateral direction. The scat- 
tering region then consists of a periodic array of paral- 
lel wires, and the lateral supercell must be chosen large 
enough to prevent an interaction between these wires. To 
limit the computational demands the supercell is chosen 
as small as possible, without a significant loss of accuracy. 

using a 



According to the results obtained in Sec. IIVBI 
2x2 supercell is already sufficient for the Na system dis- 
cussed here. We average the conductance over the 2D 
Brillouin zone (BZ) 
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where iV|| is the number of k|| -points used for the BZ sam- 
pling. Calculating the conductance for an infinitely large 
supercell would include contributions from off-diagonal 
transmission amplitudes between different k|| . From our 
results we conclude that their contribution is small as 
compared to the contribution of the diagonal terms Gk,, 
already for a 2 x 2 supercell. 

The calculated conductance Gk,, as a function of k|| for 
four and five atomic wires is shown in Fig. 1181 The dis- 
persion of the conductance is relatively small around the 
Fermi energy, which means that using a coarse k|| grid to 
calculate the conductance in the linear response regime is 
reasonable. The results discussed in sections have been 
obtained using a 6 x 6 grid in the BZ (6k||-points in 
the irreducible BZ). Fig. 1181 shows that even sampling 



the BZ with a single k|| -point can give a reasonable re- 
sult. This can be accidental, however, since the figure 
also demonstrates that the dispersion is quite large both 
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FIG. 18: (Color online) Conductance (in units of Go) as a 
function of energy for four (top) and five (bottom) atomic 
wires. The dots correspond to the conductance for various 
k||-points of a 20 x 20 grid in the BZ (55k||-points in the 
irreducible BZ). The thick line is the averaged conductance. 
The red (light grey) and blue (dark grey) dots correspond to 
k|| = (0,0) and (0.5,0.5), respectively. 



for energies lower and for energies higher than the Fermi 
energy. Especially for higher energies a single kii -point is 
clearly insufficient for calculating the conductance. This 
regime becomes important if current- voltage characteris- 
tics are calculated, because such calculations require an 
integration over a wide energy range. The large disper- 
sion at higher energies is related to an increased number 
of van Hove singularities in the leads. This suggests that 
k|| -point sampling can be important for leads containing 
atoms with a valency higher than 1 , because the number 
of van Hove singularities is then usually also higher due 
to a more complicated band structured 
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